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Fhis paper provides a prescription for the turbulent viscosity in rotating shear flows for use e.g. in 
geophysical and astrophysical contexts. Fhis prescription is the result of the detailed analysis of the 
experimental data obtained in several studies of the transition to turbulence and turbulent transport 
in Faylor-Couette flow. We first introduce a new set of control parameters, based on dynamical 
rather than geometrical considerations, so that the analysis applies more naturally to rotating shear 
flows in general and not only to Faylor-Couette flow. We then investigate the transition thresholds 
in the supercritical and the subcritical regime in order to extract their general dependencies on the 
control parameters. Fhe inspection of the mean profiles provides us with some general hints on 
the mean to laminar shear ratio. Fhen the examination of the torque data allows us to propose a 
decomposition of the torque dependence on the control parameters in two terms, one completely 
given by measurements in the case where the outer cylinder is at rest, the other one being a universal 
function provided here from experimental fits. As a result, we obtain a general expression for the 
turbulent viscosity and compare it to existing prescription in the literature. Finally, throughout all 
the paper we discuss the influence of additional effects such as stratification or magnetic fields. 



PACS numbers: 

I. INTRODUCTION 

One of the basic principle of fluid mechanics is the 
so-called "Reynolds similarity principle": no matter 
their composition, size, nature, different flow obeying 
the same equations with the same control parameters 
will follow the same dynamics. This principle has been 
used a lot in engineering to built e.g. prototypes of 
bridges to be tested in wind tunnels before construc- 
tion. To obtain easy-to-use prototypes with realistic 
control parameters, one then decreases the size but 
increases the velocity of the in-flowing wind so as to 
keep constant the Reynolds number, controlling the 
dynamics of the flow. This principle could also be of 
great interest for certain astrophysical flows, whose 
dynamics could well be approached by simple labora- 
tory flows. A good example is circum-stellar disk. In 
[l[ , it has been shown that under simple, but founded 
approximations, their equation of motions were sim- 
ilar to the equation of motion of an incompressible 
rotating shear flow, with penetrable boundary con- 
ditions and cylindrical geometry. This kind of flow 
can be achieved in the Couette- Taylor flow, a fluid 
layer sheared between two coaxial cylinders rotating 
at different speed, while penetrable boundary condi- 
tions can be obtained using porous material. On more 
general grounds, the Taylor-Couette device is also an 



excellent prototype to study transport properties of 
most astrophysical or geophysical rotating shear flows: 
depending on the rotation speed of each cylinder, one 
can obtain various flow regimes with increasing or de- 
creasing angular velocity and/or angular momentum. 

The Taylor-Couette flow is a classical example of 
simple system with complex and rich stability prop- 
erties, as well as prototype of anisotropic, inhomoge- 
neous turbulence. It has therefore motivated a great 
amount of laboratory experiments, and is even the 
topic of a major international conference. Tagg (see 
http://carbon.cudenver.edu/rtagg) has conducted a 
bibliography on Taylor-Couette flow, which gives a 
good idea of the prototype status of this flow. 

Here, we make use of the many results obtained so 
far for the Taylor-Couette experiment, regarding tran- 
sition to turbulence, or turbulence properties to pro- 
pose a practical prescription for the turbulent viscos- 
ity as a function of the radial position and the control 
parameters. It reads 



1 4 Gi{Re,T]) Si aT , 



(1) 



where Re, Rq, Rc are the control parameters (Section 
III Bp . Gi(Re,rf) is the torque measured when only 
the inner cylinder is rotating (Section [VD[1 . h(R^,r]) 
is a universal function provided in Section IV Dl and 
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Siam/ S is the ratio of the laminar to the mean shear, 
which encodes all the radial dependence as illustrated 
in Section IV Dl S and r are typical shear and radius 
of the considered flow. Most of the results we use 
here have been published elsewhere, except recent ex- 
perimental results obtained by Richard [2] . Our work 
therefore completes and generalizes the approach pio- 
neered by Zeldovich [3] . with subsequent contributions 
by [IHfi! , in which usually only one aspect of the ex- 
periments has been considered. 

An application of these findings to circumstellar 
disks using the Reynolds similarity principle can be 
found in Hersant et al. [lj thereby providing a physical 
explanation of several observable indicators of turbu- 
lent transport. 



II. TAYLOR-COUETTE FLOW 
A. Stationary flow 

The Taylor- Couette flow is obtained in the gap d be- 
tween two coaxial rotating cylinders of radii rj, ) rotat- 
ing at independent velocities fii. Q . For the purpose of 
generality and to allow further comparison with astro- 
physical flows, the velocity field at the inner cylinder 
boundary may have a non-zero radial component. 

The hydrodynamic equation of motions for an in- 
compressible flow are given by: 

1. 



dtu + u-Vu = — Vp + vAu 



Vu = 0. 



(2) 



where p and v are respectively the fluid density and 
kinematic viscosity, u is the velocity, and p is the pres- 
sure. Equation @ admits a simple basic stationary 
solution, with axial and translation symmetry along 
the cylinders rotation axis (the velocity only depends 
on r) . It is given by a flow with zero vertical velocity, 
and radial and azimuthal velocity given by: 
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where A and B are constants and a = Kjv. This 
basic laminar state depends on three constants A, B 
and K, which can be related to the rotation velocities 
at the inner and outer boundaries Ml : 
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where ry = Ti/r and a = Kjv = u r (r i )r i /i' is the 
radial Reynolds number, based on the radial velocity 
through the wall of the inner cylinder. 
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FIG. 1: Influence of the radial circulation onto the az- 
imuthal profile. Line: case a = 0; □: a = —1; ■ a = 1; 
o: a = — 10; •: a = 10. The upper panel is with 
fli/^lo = 0.86; the lower panel is with fli/£l = 0. The 
radius ratio r;/r has been arbitrarily fixed at 0.7 



The radial circulation is quantified by the value of 
a. It is positive for outward motions. For imperme- 
able cylinders, a — K = and one has the "classical" 
Taylor-Couette flows. For a porous internal cylinder, 
one obtains a Taylor-Couette flow with radial circu- 
lation. The strength of the radial circulation can be 
controlled by using more or less porous cylinders @). 
Fig. [T] provides an example of the influence of the ra- 
dial circulation on the azimuthal profile. 

In practice, even for impermeable cylinders, the flow 
is not purely azimuthal. Because of the finite ver- 
tical extent of the apparatus, a large-scale - Ekman 
- circulation is established through the effect of the 
top and bottom boundaries. This circulation depends 
on the ratio of radii and velocities, and on the top 
and bottom boundary conditions @, 0] . Its signature 
is easy to detect by profile monitoring, or by mea- 
suring the difference between the torque at the inner 
and outer cylinder [lo| . Of course this circulation is 
both radial and vertical and it varies along the cylin- 
ders axis. Also its intensity is not easy to control, 
since it is not fixed externally, but results from a non- 
trivial equilibrium within the flow. Still, at a given 
axial position, one may estimate this intensity by a fit 
of laminar profile using ([3]) and (@|. To simplify the 
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exploration of the parameter space, we shall restrict 
ourselves to the case of a ~ 0, and study separately 
the influence of this parameter. In the laboratory 
minimizing circulation effect is achieved by working 
with tall cylinders and consider only a fraction of the 
flow located at a distance to the top of about 1/3 of 
the total height, where the radial velocity is expected 
to be the weakest. Specific influence of a on stability 
and transport properties will be considered in Section 
HnDl and lVE2] 



B. Control parameters 

Dimensional considerations show that there are 
only four independent non-dimensional numbers to 
characterize the system, which can be chosen in vari- 
ous ways. 

1. Traditional choice 

The traditional choice is to consider d = r — n as 
the unit length, and d 2 jv as the unit time. With this 
choice, the dimensionless equations of motions are: 



d t u* + u*Vu* = 
V • u* = 0, 

with boundary conditions : 



-Vp* + Au* 



u*(r0 
u*(r Q ) 



where : 
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In the following, we will omit the star superscript indi- 
cating non-dimensional quantity. The present choice 
of unit amounts to define the control parameters by 
non-dimensional boundary conditions. When compar- 
ing flows that do not share identical geometry, it is of 
interest to identify control parameters characterizing 
the dynamical properties of the flows. 



2. Dynamics motivated choice 

In the case of rotating shear flows, it is convenient 
to write the equations in an arbitrary rotating frame 
with angular velocity Q r f, choose d as unit length, the 



inverse of a typical shear S as unit time and f as a 
typical radius. Furthermore, it is useful to introduce 
the " advection shear term" proposed in 0, [ll[ 

w.'V'w = w .V (w r e r ) + [rw .V (w^/r^e^+w .V (w z e z ) . 

so that the contribution of the mean flow derivative to 
the modified advection term vanishes when the flow 
is not sheared, for azimuthal axisymmetric flow. As a 
result, one has 

9 t w + w-V'w = — V7r — i?fje z x w 

+Rc ( %e r - ^e q 
\r/r r/r 



-Re' 1 Aw, 



V-w = 0, 
with boundary conditions : 



, , Rn r 

w(ri, ) = u(r i)0 ) - — re^ 
ZHc r 



(9) 



(10) 



where : 
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Rc = d/r 



(11) 
(12) 



are the dynamical control parameters for a given ra- 
dial circulation a. Re is an azimuthal Reynolds num- 
ber, measuring the influence of shear. i?n is a rotation 
number, measuring the influence of rotation. Note 
that 7r now also includes the centrifugal force term. 

In this general formulation, one is free to choose 
Cl r f. It is convenient to choose Q r f as a typical rate of 
rotation (l so that one can easily compare the Taylor- 
Couette case to the case of a plane shear in a ro- 
tating frame. For instance one can choose f2 r / so 
that wg(ri) = —wg(r ) in order to restore the symme- 
try between the two walls boundary conditions. This 
choice of fi amounts to fix f by = Q. For con- 
sistency, it is then convenient to choose S — S lam (f). 
In this context and with a = 0, it is easy to relate 
the above control parameters (Re, Rq, Rc) to the tra- 
ditional choice (Ri,R ,r]) : 



Re = — = -^—\r,R -R i \ 
v I + 77 

2fi . Ri + R 

S -qRo - Ri 



Rc = 



i-n 



1/2 ■ 



(13) 



The above control parameters have been introduced 
so that their definition apply to rotating shear flows 
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in general and not only to the Taylor-Couette geom- 
etry. It is very easy in this formulation to relate the 
Taylor-Couette flow to the plane Couette flow with 
rotation, by simply considering the limit Rq — > 0. 
Also, in the astrophysical context, one often consid- 
ers asymptotic angular velocity profiles of the form 
f2(r) ~ r~ q where q then fully characterizes the flow. 
In that case q = —dlnft/dlnr = —2/Rq, which is 
a simple relation to situate astrophysical profiles in 
the control parameters space of the Taylor-Couette 
flows. From the hydrodynamic viewpoint, an impor- 
tant characteristic of the flow profile is the sign of the 
shear compared to the sign of the angular velocity, 
which defines cyclonic and anticy Ionic flows. For the 
co-rotating laminar Taylor-Couette flow, the sign of 
the local ratio ft(r)/S(r) is constant across the whole 
flow and is thus simply given by the sign of the rota- 
tion number (Rq > for cyclonic flows and Rq < 
for anticylonic flows). Finally let us recall that an 
analogy exist between Taylor-Couette and Rayleigh- 



Benard convection (see 12[ for details and [13j for a 



review), which calls for an even larger generalization 
of the control parameters definition. 

Figure [2] displays the characteristic values taken 
by the new parameters in the usual parameter space 
(Ri,R ) for co-rotating cylinders. It also helps to sit- 
uate cyclonic and anti-cyclonic flows, as well as pro- 
totypes of astrophysical flows. 



III. STABILITY PROPERTIES 

A. Inviscid limit and data sources for viscid 
flows 

As usual when considering stability properties, one 
must distinguish stability against infinitesimal distur- 
bances - linear stability - from that against finite am- 
plitude ones - non-linear stability. When the basic 
flow is unstable against finite amplitude disturbance, 
but linearly stable, it is called subcritical, by contrast 
with the supercritical case for which the first possi- 
ble destabilization is linear (see [H], EEI for further 
details). 

In the inviscid limit (Re — » oo), and for axisymmet- 
ric disturbances, the linear stability properties of the 
flow are governed by the Rayleigh criterion. The fluid 
is stable if the Rayleigh discriminant is everywhere 
positive: 



-d r L(r) > 0, 
r 



(14) 



where L(r) 



2 il(r) is the specific angular momen- 



tum. Applying this criterion to the laminar profile 
leads to 



(R Q + l)(R u + l-f 2 /r 2 )>0. 



(15) 



Since f/r varies between 1/r] and r], one obtains that in 
the inviscid limit, the flow is unstable against infinites- 
imal axisymmetric disturbances when Rq~ < Rq < 
where = -1, respectively R^ + = I/tj-1, 

are the marginal stability thresholds in the inviscid 
limit (superscript oo) in the cyclonic case (Rq > 0, 
subscript +), respectively anticyclonic case (Rq < 0, 
subscript — ). These Rayleigh limits are also displayed 
on figure [2] where they have to be seen as asymptotic. 
As a matter of fact, this information is rather poor : 

• non-axisymmetric disturbances can be more 
destabilizing than axisymetric ones, so that the 
flow could be linearly unstable in part of the 
linearly stable domain; 

• viscous damping will probably reduce the lin- 
early unstable domain; 

• finally, finite amplitude disturbances may seri- 
ously reduce the stable domain. 

In the following, assuming that the axisymmetric 
disturbances are indeed the most dangerous one at the 
linear level - which up to now is validated both exper- 
imentally and numerically -, we will consider the two 
last items. On one side, we will review the existing re- 
sults on the effect of viscosity in the supercritical case, 
which will provide us with a critical Reynolds num- 
ber as a function of the other parameters R c (Rq,t]). 
On the other side, we will investigate the subcritical 
stability limit, when the flow is linearly stable and 
try to figure out what is the behavior of the min- 
imal Reynolds number for self-sustained turbulence 

Rg(RQ,T])- 

These boundaries can be estimated via different 
tools, depending of the type of experiment and avail- 
able measurements. In numerical experiments, the 
simplest way to estimate the stability boundary in the 
linear case is through a modal decomposition and a 
monitoring of real part of the largest eigenvalue. In 
laboratory experiments, at least three different tools 
have been used: i) torque measurements; ii) flow vi- 
sualization; iii) mean velocity profile measurements. 
Torque measurements have been traditionally used in 
the past d, El]. Their advantage is their accuracy 
and their flexibility to detect other transition at larger 
Reynolds numbers. Their inconvenience is their diffi- 
culty of implementation in the case where both cylin- 
ders are rotating. Flow visualizations allow discrimi- 
nating between laminar and turbulent flows but suffer 
from the lack of quantitative information on the flow. 
Mean velocity profile measurement is a third alterna- 
tive, which allows determination of critical Reynolds 
number from deviation of velocity profiles with respect 
to laminar value, or changes of regime. This tech- 
nique is more local in nature, and requires advanced 
techniques of in-flow measurements. In the sequel, we 
shall use data from several sources, described in ta- 
ble |TJ Except for the data of Richard, all of them have 
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FIG. 2: Parameters space and some TC flow properties for 
co-rotating cylinders (77 = 0.72). Flows with positive gra- 
dient of angular momentum L — r 2 Q, but negative (resp. 
positive) gradient of angular velocity are referred to as 
Keplerian (resp. stellar). The shaded area corresponds to 
Rayleigh unstable flows (supercritical case). 

been published. Those by Richard are available in his 
thesis manuscript [|[ . We take the opportunity of this 
synthesis to integrate them in a larger perspective. 



77 Rn source 

1 [0,0.1] Tillmark et al. [17] 

0.983 ~ Prigent et al. [18] 

0.724 -0.276 Lewis et al. [19] 

0.7 [0, 0.6] Richard [2] 

0.7 [-1.5,-1] Richard [2] 

[0.79, 0.97] I/77 - 1 Taylor [16] 



0.68, 0.85, 0.93 [-0.7, 0.5] Wendt [9] 
TABLE I: Experimental data and sources 



B. Super-critical case 

Numerous experimental set ups were used to study 
the stability boundary in the linear case, starting from 
the early experiments of Couette [2(| , Taylor [21| , and 
Donnelly and Fultz [22j. The viscosity damps the in- 
stability until Re > R c (Rn,rf), corresponding to the 
transition from the laminar flow to the so-called Tay- 
lor vortices flow. Figure [3] displays the numerical data 
by Snyder [23[ , providing the stability threshold R c as 
a function of Rq for three gap size (77 = 0.935, 0.8, 0.2), 
and illustrates the influence of the curvature on the in- 
stability threshold. The experimental data of Prigent 
et al. 18], at 77 = 0.983 is also reported. 



-1.5 -1 -0.5 0.5 1 

FIG. 3: The linear stability boundary. Numerical data by 
Snyder EB: 77 = 0.2; A: 77 = 0.8, 77 = 0.935. Experi- 
mental data by Prigent et al. V: 77 = 0.983. Continuous 
line: Lezius and Johnston plane Couette or small gap limit 
stability criteria. Dashed line : Esser and Grossman pre- 
diction for 77 = 0.8 and 77 = 0.2. 



As 77 — y 1 (rotating plane Couette limit), the sta- 
bility curve becomes symmetric around -Rq = —1/2 
and diverges at Rq = or —1. This is in agreement 
with the linear stability criterion for the rotating plane 
Couette flow [24B26J , a generalization of the first ex- 
act result giving the linear stability of the non-rotating 
plane Couette flow for all Reynolds number E3. The 
observed symmetry actually reflects symmetry in the 
rotating plane Couette. The linearized equations of 
motions are invariant by the transformation exchang- 
ing streamwise and normal to the walls coordinates 
and velocities (corresponding to exchanging r with <f> 
and u r and u$, in Taylor-Couette). This transforma- 
tion changes Rq into —1 — Rq, hence the symmetry 
around Rq = —1/2. When 77 becomes smaller than 1, 
curvature enters into play and breaks the symmetry 
resulting in less and less symmetrical curves, as can 
be observed for 77 = 0.2. 

The above stability boundary can be recovered nu- 
merically by classical stability analysis, using e.g. nor- 
mal mode analysis with numerical solutions [28]. In- 
terestingly a very good approximate analytical for- 
mula in the whole parameter space has recently been 
derived by Esser and Grossmann [29(. It is: 

R\ {Rq + l)(Rn + l- -X^J = 

with 



G 



dn _ 

d 

a(rj) 




(17) 



where A(x) is a function equal to x if x < 1 and equal 
to 1 if x > 1. Continuous lines on figure |3] give a good 
insight on the validity of the above formula. Let us 
underline that in the absence of the above formula, 
small and wide gap approximations were often used. 
Whereas the small gap approximation 
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-Rq(Rq + 1) 



(18) 



works rather well until rj = 0.8, the large gap one 
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1 



4 Q(Ra + l)' 



(19) 



where Q(x) is the Heavyside function, remains a very 
poor approximation even for 77 = 0.2. Note that 
the formula (j!6p defines two critical rotation num- 
ber for which the critical Reynolds number diverges: 
R%~ = -1 and R c n + such that R$ = -1 + 1/ -qx 2 . This 
number has been computed for various 0.7 < 77 < 1 
and is shown on figured! One sees that it is very well 
approximated by the formula R n + = (1 — rf)/r]. This 
remark is used in the next section. 

In this supercritical situation, the flow undergoes 
several other bifurcations following the first linear in- 
stability and turns into more and more complex pat- 
terns, eventually leading to turbulence. Interestingly, 
at much larger Reynolds number, an additional tran- 
sition have been reported [30j. One indeed observes 
a change in the torque dependence on the Reynolds 
number, which could be associated with a feature- 
less turbulence regime. Sometime called "hard tur- 
bulence", this regime is observed for Re > Rt- For 
reasons that will become clearer, we defer its discus- 
sion after study of the torque. 
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FIG. 4: The critical rotation number Rq~ as a function of 
the gap size rj. ■: Computed from the analytic formula of 
Esser and Grossman. Plain line: (1 — rj)/rj. 



on figure (3J giving R g as a function of Rq for different 
value of 77. 
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C. Sub-critical case 

In the absence of general theory for globally sub- 
critical transition to turbulence, the non-linear stabil- 
ity boundary has only been explored experimentally. 
Wendt [9| and Taylor [16[ consider the case with inner 
cylinder at rest, corresponding to R^ = 1/rj— 1 > 0, at 
various gap size, using torque measurements. A more 
recent experiment by Richard Q explores the domain 
— 1.5 < Rq < —1 and Rq > 0.5, at fixed gap size 
rj = 0.7, using flow visualizations. Finally one also 
has the measurement conducted in a rotating plane 
Couette flow (77 = 1) by Tillmark and Alfredsson [17j 
for Rq > 0. The corresponding results are reported 



FIG. 5: The non-linear stability boundary. Cyclonic flow: 
A: Taylor data with inner cylinder at rest and 0.7 < 77 < 
0.935; o: Richard data with 77 = 0.7; 0: Tilmark data for 
rotating plane Couette flow 77 = 1; dotted line: linear fit 
of Tilmark's data; plain line: linear fit of Richard's data. 
Anticyclonic flow: •: Richard data with 77 = 0.7; plain 
line: linear fit of Richard's data. 

One must be very cautious when looking at this 
naive representation of the data, especially on the cy- 
clonic side. First, the data are presented for different 
values of 77. Especially for the data of Taylor (A), 
each point is a different 77. The fact that the data 
look aligned through all values of 77 is an artifact of 
the representation as illustrated by the extrapolation 
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of the linear fit of Tillmark's data. Second, looking 
at this figure — seems to play a similar role in 
the cyclonic regime than Rq = — 1 in the anticyclonic 
case. As we have seen above, when studying the lin- 
ear stability, this is true for rj = 1 only. As discussed 
in previous section, the correct value of the marginal 
stability is approximately equal to the inviscid limit 



for the cyclonic case F£ 



He 



l/77-l. Tay- 



lor's data are actually given at this precise value of 
i?o, because Taylor performed his experiments with 
the internal cylinder at rest. This condition imposes 
Rq = 1/77—1, which coincides with the marginal sta- 
bility limit. In the following, we shall try to extract 
from this data the maximal knowledge about the de- 
pendence of R g ( + '' on Rq and 77, both in the cyclonic 
and anticyclonic case. 

All the data about the manifold R g ^ + '~\Ra,r]) are 
obtained close to its intersection with the manifold 
Rn = Rh i+ '~ ] - Therefore one first has to estimate the 
locus of the intersection between these two manifolds 



that is R, 



(+,-) 



H 7 /)' 7 ?) = f^ + '~\v): then the 
variation of R g + with Rq, close to the manifold, at 
the intersection. 

Let us first consider the cyclonic case. One can 
take benefit of Taylor's and Wendt's data to estimate 
f + (j]) as proposed by Richard and Zahn Q. The fact 
that the data are read from the original figure of Tay- 
lor and Wendt however induces a natural error bar 
in the determination of the critical Reynolds number, 
as illustrated on figure El where several estimates, ob- 
tained by different authors, are reported. 




FIG. 6: Subcritical thresholds in the cyclonic regime, ob- 
tained with the inner cylinder at rest, that is Rn = Rn c+ . 
□ : Wendt's data; A: Taylor's data; o: Richard's data; 
plain and dotted line: fit of f + {rj) = R g + (Rq + (rj) , rj) 
(see text for details). The size of the symbol denotes dif- 
ferent estimate by Richard and Zahn (small), Zeldovich 
(medium) and present authors (large) based on published 
figures of Taylor and Wendt. 



regime in 1 — r\ given by /+ (77) = 1400 + 550000(1 — rj) 2 
and proposed by Richard and Zahn provides a good 
upper estimate of the function. A linear trend in 1 — 77, 
with slope 136000 gives a good lower estimate of the 
data for 1 — 77 < 0.1, as shown on figure |6l Clearly, 
more precise estimate of this function using modern 
data will be welcome. Note that at 77 — > 1, the func- 
tion tends to a constant / + (1) = 1400 that is nothing 
but Rg — 1400, the global stability threshold mea- 
sured independently by Tillmark and Alfredsson [TtJ 
and Dauchot and Daviaud [3l| in the non-rotating 
plane Couette flow. The second step is to propose 
a linear development in Rq 
estimate: 



Rq + , close to the above 



f+W+a+WtRn-R^iv)). (20) 



For 77 = 1 one recovers the linear fit proposed by 
Tillmark and Alfredson (plotted and extrapolated on 
fig [5]) for the rotating plane Couette flow: 



R+(Rn,r) = l) = 1400 



26000i? o . 



(21) 



that is a+(?7 = 1) = 26000. For r/ = 0.7, the linear fit 
of Richard's data (plotted and extrapolated on fig [5]) 



leads to a+(?7 = 0.7) = 59000. 

In the anticyclonic case, the situation is simpler be- 
cause i?£ 2 ~ = — 1, does not depend on 77. On the other 
hand, data are available for a unique value 77 = 0.7, so 
that one cannot estimate f~{rj). The only fit that can 
be performed in this state of experimental knowledge 
is: 



R7(R Q ,0.7) = /-(0.7)+a"(0.7) | Rq-R? 



(22) 



One finds /"(0.7) = 1300 and a" (0.7) = 21000 and 
the fit is displayed on figure [5j In the anticyclonic 
regime, at least for this value of 77, one recovers a 
dependence on the rotation similar to that of the plane 
Couette flow. Also remarkable is the fact that /~ is 
so close to Rg C in the non-rotating case. 

Altogether the data collected to date suggest that, 
in the linearly stable regime, the Reynolds number 
of transition to subcritical turbulence be well repre- 
sented by 



f ± (ri) + a ± (r ] )\R n -R 



C,±l 



(23) 



with 1 x 10 5 (1 - 77) 1 < f+{rf) < 1400 + 5.5 x 10 5 (1 - 
?7) 2 ,/~(0.7) = 1300 and 21000 < a ± < 59000. It is 
difficult to distinguish the effect of experimental pro- 
cedures from the effects of gap width dependence in 
the present parameter range. 



D. Influence of radial circulation 



Super- critical case 



Because of this error, it is difficult to give a precise 
fit of the function /+. One sees that the quadratic 



The influence of radial circulation on the linear sta- 
bility onset has been studied numerically by Min and 
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Lueptow [8J. They observed that an inward radial 
flow and strong outward flow have a stabilizing effect, 
while a weak outward flow has a destabilizing effect. 
We may use their data to get more precise estimates 
for the case a = —3/2, (q = 3/2, Keplerian case). Fig- 
ure [7] shows the ratio R c (a = — 3/2) /R c (a = 0) — 1 
as a function of rj for f2 Q /f2j = 0. One sees that the 
variation is quasi-linear. 
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(b) 
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FIG. 7: R c (a = -3/2)/R c (a = 0) - 1 as a function of (a) 
77 for fio/ft = 0; (b) n /fij = for 77 = 0.85. □ : data 
from Min and Lueptow. The dotted lines are the fit eq. 
423 and J25). 



2. Sub-critical case 

The influence of the radial circulation on the non- 
linear stability has not been systematically studied. 
However, we can get partial answers from the exper- 
iments of Wendt and Richard 0], where the in- 
fluence of the top and bottom circulation on the on- 
set of stability has been studied. Both Richard and 
Wendt investigated the stability boundary with dif- 
ferent boundary conditions. One boundary condition 
was with the bottom attached to the outer cylinder. 
In this case, the circulation is mainly in the anti- 
clockwise direction, with radial velocities outwards at 
the bottom (a > 0). Another boundary condition 
was with the bottom attached to the inner cylinder 
(at rest). In that case, the circulation is in the op- 
posite direction, with inward radial velocities at the 
bottom (a < 0). A last boundary condition was in- 
termediate between the two, with only part of the 
bottom attached to the outer cylinder. In neither 
case, noticeable change of the stability boundary has 
been noticed, which means that at this aspect ratio 
7; = 0.7, the radial circulation induced by the bound- 
ary conditions has an impact on the subcritical thresh- 
old Reynolds number which is less than 10 per cent 
(accuracy of the measurements). 



E. Influence of aspect ratio 

Most of the experimental set-ups described in this 
paper have a very large aspect ratio V = H/d ^> 1. 
Keplerian disks are characterized by a small aspect 
ratio H/d ~ 0.01 — 0.1. It would be interesting to 
conduct systematic studies of the variation of T onto 
the stability and transport properties. The influence 
of r onto the instability threshold, in the case of 
outer cylinder at rest has been computed by Chan- 
drasekhar [H|, Snyder This is illustrated in Fig. 
[H The critical Reynolds number is increased, as T is 
decreased. It follows an approximate law: 



fl c (r) = R c (1 + r 



0(T~ 



(26) 



A best fit gives: 

^7-^ = 1+0.12(1-,), £ = 0. (24) 

R c (a = 0) Qi 



On the same graph, we show R c (a = — 3/2) / Re c (a — 
0) — 1 as a function of f2 /^i for 77 = 0.85. A best fit 
gives: 



Re + (a = -3/2) 
Re + (a = 0) =1 + 1 - 16 U' ' " 



(25) 



This behavior can be understood if one says that as 
r becomes smaller, the smallest relevant length scale 
in the problem become H instead of d. The relevant 
Reynolds number has thus to be corrected by a factor 
(H/d) 2 , hence, the T -2 law. However, another ex- 
perimental study by Park et al. [32J suggests that the 
physical relevant length scale is V Hd instead of H. A 
possible explanation of the difference is through the 
Ekman circulation, which is present in experiments 
and not in numeric. This circulation may couple verti- 
cal and radial velocities, leading to an effective length 
scale. The only way to settle this issue is through 
smaller aspect ratio systematic laboratory and numer- 
ical experiments. 
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FIG. 8: Influence of aspect ratio onto the critical Reynolds 
number for instability, in the case where the outer cylinder 
is at rest. □: numerical data by Chandrasekhar. The 
dotted line is a power law fit 0.9r -2 . 

F. Structural stability 

To close this section, it is interesting to consider 
the influence of additional physical forces that may 
be relevant to astrophysical flows. In the sequel, we 
only give a summary of the main experimental or the- 
oretical results obtained, referring to the publications 
for more details. 



1. Magnetic field 

The influence of a vertical magnetic field on the 
stability of a Taylor-Couette flow has been studied 
theoretically |28l [33J and experimentally by Donnelly 
and Ozima |34| using mercury. Applications to as- 
trophysics have been discussed by Balbus and Haw- 
ley 35]. This motivated a lot of numerical work on 
this instability. For references, see e.g. [36j ■ 

In the inviscid limit, the presence of a magnetic 
field changes the Rayleigh criteria (jT4)) . For exam- 
ple, in the case of a magnetic field given by B = 
\/Mp(0i Hfi(r), Hz), the sufficient condition for sta- 
bility is now [38l . |39j : 

r 2 8 r n 2 - ^d r {r 2 Hl) > 0. (27) 

Therefore, anti-cyclonic flow, with i?o < are now 
potentially linearly unstable in the presence of a 
magnetic field with no azimuthal and radial compo- 
nent |35j |. 

The linear instability in the presence of dissipation 
has only been studied numerically. A first observa- 
tion was that boundary conditions (e.g. insulating 
or conducting walls) are relevant to determining the 
asymptotic behaviors [28j |. The proposed explanation 
is that the magnetic field makes the flow adjoin the 



walls for longer distances, so that the viscous dissi- 
pation remains comparable to the Joule dissipation 
at all fields. A second observation is the importance 
of the magnetic Prandtl number P m = v / ' K m (n m is 
the magnetic diffusivity) on the instability [361 |3jj . 
On general grounds, it seems that at small Prandtl 
numbers, the magnetic field stabilizes the flow in the 
supercritical case, while at large Prandtl numbers, the 
magnetic field destabilizes the flow. In the subcritical 
case, the magnetic field can excite a linear instability 
for anti-cyclonic flow, at any Prandtl number. This is 
illustrated in Fig. [5] 

Scaling of critical Reynolds number with magnetic 
Prandtl numbers have been found: in the supercrit- 
ical case — 1 < Rn < 0, the critical Reynolds num- 
ber scales like Pm 1 ^ 2 [37[. In the subcritical case 
Rq = —I, the critical Reynolds number scales like 
Pm- 1 . 
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FIG. 9: Influence of body forces on the linear stability 
boundary. EB: without forces, r\ — 0.2, numerical data 
from Snyder (1968). With vertical constant magnetic field, 
at Pm = 1 (□) and Pm = 1CT 5 (■), 77 = 0.27; numerical 
data from Riidiger et al, 2003; V: with vertical stratifica- 
tion, r\ = 0.2; data from Whithjack and Chen (1974). 



2. Vertical stratification 

A vertical stable stratification added onto the flow 
plays the same role as a vertical magnetic field at low 
Pm. In the inviscid limit, its presence changes the 
Rayleigh criteria into r 2 d r tt 2 > [H 52) ■ This 
means that all anti-cyclonic flows are potentially lin- 
early unstable. The role of dissipation on the insta- 
bility has been studied numerically [U |42| ) and ex- 
perimentally [43|, IHI- It was found that stratifica- 
tion stabilizes the flow in the GSPC regime, while it 
destabilizes it in the GSBC anti-cyclonic regime. The 
critical Reynolds number was found to scale with the 
Froude number (ratio of rotation frequency to Brunt 
Vaissala frequency) like Fr~ 2 , and to scale with the 
Prandtl number (ratio of viscosity to heat diffusivity) 
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like Pr- 1 ' 2 . 



3. Radial stratification 

A radial temperature gradient applied to the flow 
changes the stability. In the inviscid limit, the 
Rayleigh criterion is modified by the radial temper- 
ature gradient into 0, |46[ : 

% ri r*n)(l-^-)-°frtf>Q, (28) 

where 1/To is the coefficient of thermal expansion and 
AT is the temperature difference between the cylin- 
ders. The last term in (|28|) induces an asymmetry be- 
tween the case with positive AT and negative AT. An 
experimental study by Snyder and Karlsson [U helps 
to quantifying the role of dissipative processes. It was 
found that both positive and negative AT have a sta- 
bilizing effect when AT is small, and a destabilizing 
effect when AT is large. A more complete exploration 
of the parameter space would be welcome, since astro- 
physical disks arc likely to be subject to this kind of 
stratification. 



4- Summary 

These studies point out an interesting dissymmetry 
between the case Rq > (cyclonic flows) and Rn < 
(anti-cyclonic flows). In many instances, the regime of 
linear instability is extended by the large scale force 
into the whole domain Rq < 0. As a result, in the 
anticyclonic regime one often has to deal with a com- 
petition between a linear destabilization mechanism 
induced by the large scale effect and the subcritical 
transition controlled by the self-sustained mechanism 
of the turbulent state. 



IV. MEAN FLOW PROFILES 

A. Supercritical case 

Turbulent mean profiles have been measured re- 
cently for different Reynolds number by Lewis and 
Swinney [l9j in the case with outer cylinder at rest. 
They observe that the mean angular momentum L = 
rue is approximately constant within the core of the 
flow: L ~ 0.hrjQ,i for Reynolds numbers between 
1.4 x 10 4 and 6 x 10 5 . At low Reynolds number, this 
feature can be explained by noting that reducing the 
angular momentum is a way to damp the linear in- 
stability, and, thus, to saturate turbulence. At larger 
Reynolds number, however, one expects the turbu- 
lence to be sustained by the shear in the same way 



as it is when there is no linear instability at all. Ac- 
cordingly, this constancy of the angular momentum 
is quite a puzzling fact. Some understanding of this 
behavior can be obtained by observing that the mean 
profiles obtained by Lewis and Swinney are actually 
in good agreement with a profile obtained by Busse 
upon maximizing turbulent transport in the limit of 
high Reynolds number [481 149|: 

This profile bears some analogy with the laminar pro- 
file, which reads: 

4 am (r) = -V 1 ^ + r(n + ls\. (30) 

In the Busse solution, the shear profile S°°(r) = 
rjf 2 S/-ir 2 = l/AS lam (r). This ratio is analog to the 
value observed at very large Reynolds number in the 
non-rotating plane Couette flow [50|. It is therefore 
a clear signature of the shear instability, with no dis- 
cernable influence of rotation, at least for the limited 
value of the rotation number (of the order of —0.28) 
considered by Lewis and Swinney. So it is interesting 
to test the Busse asymptotic profile using other data, 
with different rotation number. This will be the pur- 
pose of the next section, where Richard data will be 
used. 

We may however not conclude this section without 
noting an intriguing property of the Busse solution. 
Considering R n turb = 2ft 00 (r)/S 00 (r), we get from 
([291): 

ifo^6 =4flQ + 3 (l-^ (31) 
(1 + rf) 

So the condition -1 < R n turb < I/77 — 1 ("linear 
stability of the turbulent profile" ) is satisfied provided 
Rn follows: 

2 — 377 + 2/7 2 1 - X] (1 + T}) 2 - 377(1 - r/) 2 

2(1+ V 2 ) < n< ^l (1+7? 2 ) ' 

(32) 

that is, in the small gap limit, —1/4 < Rq < 0. As 
we shall see in the sequel, this is precisely the range 
of value where the torque is extremum. 



B. Subcritical case 

For the turbulent flow following the subcritical tran- 
sition, we use the data of Richard [2j , collected for dif- 
ferent Reynolds numbers and rotation numbers. Fig- 
ure [TU] displays typical turbulent mean profiles in both 
the cyclonic and anticyclonic cases, for comparison 
with the laminar and the Busse profiles. 



11 




0.75 0.8 0.85 0.9 0.95 1 




FIG. 10: Mean velocity profiles from Richard at Re/R g = 
1.6: (a): cyclonic case (Rq — 0.39); (b): anti-cyclonic case 
(Rn = —0.6). Dotted line: laminar profile. Continuous 
line: Busse solution eg. (|29[) . 



One notices the profile tendency to evolve from the 
laminar one to the Busse solution, even if they are still 
very far away from the extremizing solution. In order 
to evaluate how fast the convergence occurs, figure ITT1 
displays the ratio of the turbulent mean shear to the 
laminar shear, both estimated at f, i.e. S(f)/S, as a 
function of the ratio of the Reynolds number to the 
threshold for shear sustained turbulence i.e. Re/R g . 
One may indeed observe a tendency of shear reduction 
as the Reynolds number increases, with a more rapid 
reduction for rotation number closer to 0. However, 
none of the case studied by Richard approaches the 
value 0.25 predicted by Busse. It would be interesting 
to conduct higher Reynolds number experiments at 
large value of the rotation number, to check whether 
rotation merely slow down the convergence towards 
the 0.25 value, or change it into a number depending 
on the rotation number. 

Also one may notice that the decrease of S(r)/S 
with Re/Rg is much faster for cyclonic flows than for 
anticyclonic ones. Figure [T2l may provide some hints 
on the origin of this dissymmetry. The first one is 
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FIG. 11: Ratio of the turbulent mean shear to the lami- 
nar shear variation with the Reynolds number (data from 
Richard). Cyclonic case: 0: Rn = 0.39; o: Rq £ 
[0.64,0.52]. Anticyclonic case: •: Rn G [-2,-1]. 



obtained by studying the radial variation of the ra- 
tio S lam /S at a given Re/R g , for different rotation 
number. This quantity provides the radial variation 
of the turbulent viscosity and thus is a good tracer 
of transport properties. One may observe an inter- 
esting tendency for cyclonic flow to display enhanced 
(resp. depleted) transport at the inner (resp. outer) 
core boundary, while anti-cyclonic flow rather displays 
depleted transport at the center, and enhanced trans- 
port at both boundaries. 

The second one is provided by the function: 



q(r) 



2fi(r) _ din ft 



dlnr 



(33) 



which may be viewed either as a local mean angular 
velocity exponent, or a local mean rotation number. 
This local exponent also plotted on figure [121 for dif- 
ferent rotation number, at Re/R g = 1.6. One clearly 
observes a tendency towards constancy of this local 
exponent in the core of the flow and a bimodal be- 
havior: cyclonic flow scatters towards q = 0.5 while 
anti-cyclonic flow scatters towards q = —1.5. We 
have observed a persistence of this behavior at larger 
Reynolds number (up to at least Re/R g = 20). 



TORQUE MEASUREMENTS AND 
TRANSPORT PROPERTIES 



The turbulent transport can be estimated via the 
torque T applied by the fluid to the rotating cylinders. 
Traditionally, one works with the non-dimensional 
torque G — Tj phi/ 2 [3(J. For laminar flows, one can 
compute this torque analytically using the laminar 
velocity profile. It varies linearly with the Reynolds 
number. 



G la 



2- 



(1-^ 



-r\Re. 



(34) 
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FIG. 12: Mean profile at Re/R g = 1.6 for (a) turbulent 
transport, traced by the ratio S larn /S; (b) local rotation 
number. T: Rn = -1.31; •: Rq = -1.41; o: R n = 0.39; 
0: Rn = 0.51. Data are from Richard. 



thing with exponent closer to 1. This regime has been 
observed up to the highest Reynolds number achieved 
in the experiment (of the order of 10 6 ). 

The experiment with inner cylinder rotating only 
covers flows such that Rq = T) — I. To check whether 
this kind of measurement is typical of torque behav- 
iors in the globally supercritical case, one must rely on 
experiments in which the outer cylinder is also in ro- 
tation. Unfortunately, the only torque measurements 
available in this case are quite older [9( and not as de- 
tailed as in the case with inner cylinder rotating. Most 
specifically, they do not extend all the way down to 
the transition region between laminar and turbulent. 
In several instances in which large Reynolds number 
are achieved, however, one may observe a steepen- 
ing of the relative torque towards the G/G lam ~ Re 
already observed in the case with inner cylinder ro- 
tating. On other measurements performed at lower 
Reynolds numbers, the relative torque displays a be- 
havior more closely related to the intermediate regime, 
with G/G lam ~ Re 1 / 2 . Altogether, this is an indica- 
tion that in the globally supercritical case, the torque 
follows three regimes: 



G 
G 
G 



aRe, Re < R r 



/? su Pi?e 3/2 , R c <Re<R Tl 
up Re 
2ttt] 



7 sup i?e 2 , Re>R Tl 



(35) 



(1-r?) 2 ' 

where (3 sup and j sup are constants to be specified later 



When the turbulence sets in, the torque applied to the 
cylinders tends to increase with respect to the laminar 
case. A good indicator of the turbulent transport can 
then be obtained by measuring G/G lam . 



A. Super-critical case 

As noticed by Richard and Zahn ||, most of the 
torque measurements available in the literature con- 
cern the case with the outer cylinder at rest (see 
e.g. [Ill H(| and references therein). In that case, we 
note that \Rq\ = \rj— 1| < 1. An example of the varia- 
tion oiG/G lam with Reynolds number is given in Fig- 
ure !13[ in an apparatus with rj — 0.724. One observes 
three types of behaviors: below a Reynolds number 
R c , i.e. in the laminar regime G/G lam = 1. Above R c , 
one observes a first regime in which G/G lam varies ap- 
proximately like a power-law, with exponent 1/2. In 
this regime, Taylor vortices can often be noticed. This 
regime continues until Re = Rt, where the torque be- 
comes stronger, and the power-law steepens into some- 



B. Sub-critical case 

The only measurements of torque in the subcritical 
case were performed by Wendt Q and Taylor [l6[ in 
experiments with the resting inner cylinder, and rotat- 
ing outer cylinder. Wendt 's experiments cover three 
different values of rj, Taylor's cover eleven values of 
rj. Taylor measurements cover sufficiently small value 
of Reynolds number so that one can see that above 
a critical Reynolds number Re = R g , the torque bi- 
furcates from the laminar value towards a regime in 
which the relative torque G/G lam behaves like Re. An 
example is given in figure [TBI Measurements by Wendt 
at larger Reynolds number display no evidence for an 
additional bifurcation. So, in the subcritical case, the 
torque presumably follows only two regimes: 

G ~ aRe, Re < R g , 

G ~ Y ub Re 2 , Re>R g , (36) 

where r ) sub is a constant that we specify in the next 
subsection. 
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FIG. 13: Relative torque G/G as a function of the 
Reynolds number. T: super-critical case with outer cylin- 
der at rest; Rn = —0.276; 77 = 0.724. Data are from Lewis 
and Swinney. □: sub-critical case, with inner cylinder at 
rest. Rn = 0.47; 77 = 0.68. Data are from Wendt. 



C. Connecting torque and thresholds 

Invoicing the continuity of the torque as a function 
of the Reynolds number at the transitions allows to 
determine the prefactors j3 sup , ^ sup and 7 s uh . In the 
supercritical case, one obtains: 



= aR- 1 / 2 , 



3 sup R 



-1/2 



\J R c Rt 



and in the subcritical case: 
j sub = 



R„ 



(37) 



(38) 



where a is known through (1341) . This enables the 
knowledge of the torque as a function of R c and Rt, or 
R g which then encode all the dependencies on Rn and 
77. This would be of great practical interest, and a pos- 
teriori gives all its importance to the work conducted 
in section ITTll since torque measurements are usually 
more difficult to perform than thresholds estimations, 
especially when both cylinders are rotating. Our ar- 
gument is admittedly very crude, so it is important to 
test its validity on available data. Figure [TJ] shows the 
comparison between the real non-dimensional torque 
measured in experiments, and the torque computed 
using only the critical Reynolds number. At low 
Reynolds number, there is a fairly large discrepancy 
but at large Reynolds, the approximate formula pro- 
vides a good estimate. 

Comparing ([37]) and (|38|) suggest to introduce 
R g sup = \/R c Rt- This new threshold, defined in the 
supercritical case, would correspond to the Reynolds 
number above which turbulence is sustained by the 
shear mechanism, and not anymore by the linear in- 
stability mechanisms. A physical basis for this ex- 
pression could be given using the observation that the 



FIG. 14: Relative torque G/G as a function of the 
Reynolds number, compared with its determination using 
critical Reynolds numbers. ▼: globally super-critical case 
with outer cylinder at rest. Tin = —0.276; 77 = 0.724. 
(Lewis and Swinney). □: globally sub-critical case, with 
inner cylinder at rest. Rn — 0.47; 77 = 0.68 (Wendt 
1933). Short dashed line: ^Re/R c with R c = 90; dot- 
dashed line: Re/y/RcRr with \J R c Rt = 957; long-dashed 



line: Re/R g where 



Rg — 



32688. The critical Reynolds 



numbers have been computed using results of section IIIII 



transition occurs in a turbulent state, where transport 
properties are augmented with respect to a quiescent, 
laminar case, in which all transport is ensured by vis- 
cous processes. This results in a delayed transition 
to the ultimate state, since the viscosity is artificially 
higher by an amount Ut/v, where v t is the turbulent 
viscosity. Using v t /v — G/Gi am , we thus get from 
"?5l) and (|3T|) an estimate of the relevant threshold as: 



R s a wp = Rt — 
y v t 



R T - 



VRTRt. (39) 



At this stage of the analysis, R g , R g up and i?+ respec- 
tively define a function of Rn and 77 on the intervals 



Rn < Rr 



< Rn < R Q + and R?+ < R Q , 



where we recall here that R n 



-1 and R n + 



I/77 — 1 . Further indication of the relevance of R s g up 
is provided by the continuity of this function with Rq 
throughout the super/sub-critical boundaries. This is 
illustrated on figure 1151 where the continuity is ob- 
tained on the cyclonic side between the Tillmark's 
data (77 = 1) and Wendt's data (77 = 0.935) and on 
the anticyclonic side between Richard's data (77 = 0.7) 
and Wendt's data (77 = 0.68). 



D. How to use torque with resting outer 
cylinder 

Torque measurements described in previous section 
suggests that at large enough Reynolds number, an 
" ultimate" regime is reached with quadratic variation 
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FIG. 15: R g as a function of Ra and 77. Anticyclonic side: 
• : Richard data (77 = 0.7); ■: Wendt data (77 = 0.68). 
Cyclonic side: (): Tilmark data (7/ = 1) and fit as in section 
UTTl □: Wendt data (77 = 0.935); o: Richard data (77 = 0.7). 
The lines are guides for the eyes to underline the continuity 
across the supercritical to subcritical domains for similar 
values of 77. 



FIG. 16: Relative torque G/d as a function of Rq and 
77. m, (resp. o): estimation from Richard data (77 = 0.7), 
based on critical Reynolds numbers, computed using re- 
sults of Section IIIII in the anticyclonic (resp. cyclonic 
case); ■, (resp. □): Wendt data (77 = 0.68,0.85,0.935), 
the square size increasing with 77 in the anticyclonic (resp. 
cyclonic case). 



with Reynolds number. This suggests that in this 
regime, the interesting parameter is the ratio of the 
torque in any configuration, to the torque measured 
in a special case. Because the case with resting outer 
cylinder is the most studied, it is of practical interest 
to choose this case as the reference, so that the rele- 
vant ratio is G/Gi, where Gi is the torque when only 
the inner cylinder is rotating. Given the above sub- 
sections, G/Gi is only a function of Rq and 77 given 
by h(Rn,r]) = R g (Rn(R = 0),i])/R g (Rn,ri), where 
R g is the generalized threshold defined in the previous 
section and displayed on figure Qj)J Figure [T5] indeed 
shows the ratio G/Gi, for different values of 77, as a 
function of the rotation number. The measurements 
for —0.8 < Rq < 0.5 are direct measurements from 
the Taylor and Wendt experiments. The measure- 
ments for Rq < — 1 and Rq > 0.5 are indirect mea- 
surements, coming from the experiment by Richard, 
in which only critical numbers from stability were de- 
duced. In that case, the torques have been computed 
using the results of previous section. All these results 
show that the non-dimensional torque behaves as: 



G{Re, Rq, 77) = G l (Re, n)h(Ra,v), 



(40) 



where Gi is the torque when only the inner cylinder 
is rotating, and h(RQ,rf) is the function of figure [TBI 
This universal function is very interesting because 
it provides good insight about the influence of the 
rotation and curvature on the torque. For rotation 
number —0.2 < Rq < 0, the torques are maximal 
and equal to the torque measured when only the in- 
ner cylinder is rotating. For rotation numbers out- 
side this range, torques tend to decrease, with a sharp 



transition towards a constant of the order of 0.1 on 
the side Rq > 0. On the other side, the transi- 
tion is softer, with an approximate quadratic inverse 
variation until the smallest available rotation number 
Rq = —1.5. From a theoretical point of view, the 
asymmetry could be linked with the different stability 
properties of the flow on either side of the curve: for 
— 1 < Rq < I/77— 1, the flow is linearly unstable, while 
it becomes liable to finite amplitude instabilities out- 
side this range. The variation we observe can also be 
linked with experimental studies by Jacquin et al [5lj , 
re-analyzed by Dubrulle and Valdettaro (52j. They 
show that rotation tends to inhibit energy dissipation 
and observed simple power laws linking the energy dis- 
sipation with and without rotation as cq^o = eji=oi?0' 
where Rq is a rotation number based on local shear 
and rotation. 

Finally, the previous discussion shows that the 
knowledge of the torque in the case with resting outer 
cylinder as a function of Re and 77 is an essential data 
to compute the torque in any other configuration. A 
theoretical model of the torque in that configuration 
has been proposed by Dubrulle and Hersant [12j . in 
the case where the boundary conditions at the cylin- 
der are smooth. It gives: 



Gi 

Gi = 0.33 



(3 + 7 7 ) 1 / 4 (7 7 7?e) 3 / 2 



(l-77) 7 /4(l +7? )l/2^ 

(3 + r,) 1 / 2 



R c <Re< R T 
(7/i?e) 2 



0.0001 



(1 - 77) 3 / 2 (l + 77) (\ n [(K(r])(r]Re) 2 }f /2 
Re > Rt, 
(I-77X3 + 77) 



(1+77)2 



(41) 
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The quality of the fit can be checked on Fig. [171 For 
Re < R c , the flow is laminar and the transport is 
ensured only by the ordinary viscosity. 




smooth boundaries 



FIG. 17: Influence of boundary conditions on torque. Case 
with two rough boundaries at r\ = 0.724, A, data from Van 
den Berg et al.; at n — 0.625, EH, data from Cadot et al.. 
The continuous lines are the formula (|43p . Case with two 
smooth boundaries at r\ — 0.724, V. Data from Lewis 
and Swinney. The dotted and the dashed-dot lines are the 
formulae (|4*T1) . 



E. Towards extended Reynolds similarity 

The link between torque and critical Reynolds num- 
ber has a powerful potential for generalization of the 
torque measurements performed in the laboratory for 
astrophysical or geophysical flows. Indeed, all the ad- 
ditional complications studied so far (aspect ratio, cir- 
culation, magnetic field, stratification, wide gap limit) 
have been found to shift the critical Reynolds number 
for linear stability by a factor function of this effect, 
like R c (effect ^ 0) = R c (effect = Q)F. Depend- 
ing on the situation, F can be interpreted as either 
a change in the effective viscosity (magnetic field), 
or a change in the effective length scale (aspect ra- 
tio, wide gap). If, on the other hand, the scaling of 
the torque with Reynolds number (i.e. the shear) re- 
mains non-affected by such a process, the computation 
done in section IV Dl are easy to generalized through 
an effective Reynolds number Re e '' = Re/F. Specif- 
ically, everything that has been said for the torque, 
in the ideal Taylor-Couette experiment, will still be 
valid with additional complication provided one re- 
places the Reynolds number by an effective Reynolds 
number, taking into account the stability modification 
induced by this effect. This principle is by no mean 
trivial and must be used with caution, even though it 
may appear as nothing more that an extension of the 
Reynolds similarity principle. In fact, it has been val- 
idated so far only in the case with vertical magnetic 



field, where it has been indeed checked by Donnelly 
and Ozima [13] that the torque scaling is unchanged 
by the magnetic field. In the sequel, we shall use this 
procedure in disks, because we noticed that it gave 
the most sensible results. It would however be impor- 
tant to check experimentally this "extended Reynolds 
similarity" principle. 



1. Influence of boundary conditions 

Experimental investigation of the Taylor-Couette 
flow with different set-up has shown that boundary 
conditions have an influence on the torque. More pre- 
cisely, it has been shown that the inclusion of one [53[ 
or two [HI rough boundary condition, in configura- 
tion with outer cylinder at rest, increases the torque 
with respect to the case with two smooth boundary 
conditions, at large Reynolds numbers. In convec- 
tive flows, a similar increase of transport properties 
is observed when changing from no-slip to stress-free 
boundary conditions |55| . In both cases, the increase 
occurs so as to increase the agreement between the 
observed value, and a value based on classical Kol- 
mogorov theory. A theoretical study of Dubrulle [56] 
explains this feature through the existence or absence 
of logarithmic corrections (see formula |4T]-b)) to scal- 
ing generated by molecular viscosity and large-scale 
velocity gradient in the vicinity of the boundary. Ob- 
viously, in the presence of a rough boundary, or under 
stress-free boundary conditions, mean large-scale ve- 
locity gradients are erased near the boundary, and no 
logarithmic correction develops. 

For two rough boundary conditions, Cadot et al [54[ 
measure Gt ~ 0.22 — 0.3i?e 2 for rj = 0.625, while van 
den Berg et al. [^observe G t ~ 0.43i?e 2 for rj = 0.73. 
The analogy with thermal convection [l2j suggest that 
Gi depends on its laminar value and on rj and Re like 



G, 



G 



lam 



X /(3 + 7 ? )(l-7 ? ) 
trough (l+rj) V ^ ' 



Using (J34J), and the experimental law, we find 
Irough = 0.017, so that: 



Gi = 0.107 



(1 + 7 ? )(1- r?) 3 / 2 



(rjRe) 



(43) 



The comparison between this formula and the exper- 
iments is made in figure 1171 For reference, we also 
added the torque in the case of two smooth boundary 
conditions, as given by (TTTj-b). 

We do not have any theory for the case with asym- 
metric boundary conditions (one rough, one smooth) . 
Laboratory experiments show that the torque lies in 
between the curve for two smooth boundary condi- 
tions and the curve for two rough boundary condi- 
tions. The exact location however depends on local 



16 



conditions in a non-trivial way (for example it is dif- 
ferent when the rough conditions applies to the (rotat- 
ing) inner cylinder or to the (resting) outer cylinder). 
The present experimental evidence therefore only al- 
lows the torque measurements with two smooth (resp. 
two rough) boundary conditions to be considered as 
lower (resp. upper) bounds for the torque, in case of 
complicated boundary conditions. 



2. Influence of radial circulation 

Torque measurements by Went, for different ge- 
ometry show that the circulation can have an in- 
fluence on the transport properties. Specifically, it 
has been observed that an outward circulation tends 
to increase the torque applied on the inner cylinder, 
while an inward circulation tends to decrease this 
torque. The difference can be quite important. At 
large Reynolds numbers, the relative increase of the 
torque G(a < 0)/G(a > 0) can be computed as a 
function of 77. This is shown in Fig. [T5] One observes 
a quasi-linear variation: 



Gja < 0) 
G(a > 0) 



12.45(1 - 77), Qi = 0. 



(44) 



The case with intermediate boundary conditions (pre- 
sumably a close to zero) lays about half way in be- 
tween the two cases so that: 



G(a < 0) 
G(a = 0) 



4.75(1 - 77), fi« = 0. 



(45) 



Close to the transition threshold, there is also an 
asymmetry between the two circulation regimes: out- 
ward circulation enhances the torque with respect to 
the laminar regime, while inward circulation decreases 
this torque! This puzzling aspect has been explained 
by Coles and Van Atta [1(| ; in absence of circulation, 
in stationary state, the torque exerted at the inner 
and outer cylinder must balance. In the presence of 
circulation, the transport of fluid toward or away from 
the plane of symmetry induces an imbalance of the two 
torques, which ceases to be equal. Coles and Van Atta 
measured this imbalance as a function of the Reynolds 
number for the case with inner cylinder at rest, at 
j] = 0.89, and with boundary conditions favoring an 
outward circulation. One observes an imbalance of 
the order of 30 to 50 percent on Fig. (fl~8|). with the 
torque on the inner cylinder being larger. These obser- 
vations suggest the following model: in the presence 
of a radial circulation, the inner and outer torques are 
modified into: 

G D (a) = G(a = 0)(l-aA(Re,r))), 

Gi(a) = G(a = 0){-l-aA(Re,r))), (46) 

where A is a positive function of 77 and Re. So, when 
reversing the circulation (going from a > to a < 
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FIG. 18: Influence of radial circulation on torque, (a) 
G(a > 0)/G(a < 0) as a function of 77 for = in 

the experiment of Wendt, at large Reynolds numbers. The 
squares are the data. The line is the fit eq. ((44} (b) Ratio 
of torque applied to outer cylinder vs. torque applied to 
the inner cylinder in the case of an outward circulation, 
77 = 0.89. 



0, the torque exerted at the inner cylinder decreases 
(in absolute value), like in Wendt data. Moreover, 
these data indicate that at large Reynolds number, 
the function A becomes independent of Re. Note also 
that according to this model, we should have Gi(a < 
0)/G l (a > 0) = -G (a > 0)/G t (a > 0. At Re = 
10 5 , 77 = 0.8, the data of Wendt provide a value of 
0.25 for this ratio, in good agreement with the value 
observed by Coles and Van Atta, see Fig. [THJ In 
this model, the total torque is zero (conservation of 
total angular momentum) only when considering the 
torque applied by the circulation on the top and the 
bottom boundary. This means that in the presence of 
a radial circulation, a non-negligible torque is likely 
to apply at the vertical boundary. This observation 
may be relevant to astrophysical disks, and jet-like 
phenomena. 
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3. Influence of a vertical magnetic field 

The influence of a constant vertical magnetic field 
on the torque has been studied by Donnelly and Oz- 
ima 34]. The measurements have been performed 
in the linear instability regime, with outer cylinder 
at rest. It is observed that an increasing magnetic 
field reduces the torque, so as to conserve the Re 3 / 2 
scaling observed at zero magnetic field (section IV Aft . 
The torque reduction is thus a function only of a non- 
dimensional magnetic field, and of rj. Examples are 
provided in Fig. HH for gap sizes 0.901 and 0.995 and 
Reynolds number Re ~ 2100. 
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FIG. 19: Influence of vertical magnetic field on torque 
in the case with inner cylinder rotating. Ratio of torque 
with magnetic field to the torque without magnetic field 
in same experimental conditions, as a function of Q, the 
non-dimensional magnetic number. The symbols are the 
data. The lines are the fit (|47[). Data are form Donnelly 
and Ozima. EB: 77 = 0.901, Re = 2105; the constant for the 
fit are b = 0.803 and c = 0.001; 0: r] = 0.995, Re = 2050; 
the constants used for the fit are b = 0.920 and c = 0.002. 

The torque reduction can be quantified by the dimen- 
sionless number Q — /j, 2 H 2 ad 2 / pv 1 where \i is the 
permeability, a is the electrical conductivity and H is 
the magnetic field in tesla. It seems to follow a simple 
law: 



b(r,) 



G l (Q = 0) X /1 + C (,,)Q : 



(47) 



where b(r]) and c(rj) are functions of the gap size. 
Physically, this torque reduction may be due to the 
elongation of the cellular vortices which occurs as 
the magnetic field is increased [34j. Mathematically, 
the reduction can be understood using the connection 
between torque and critical number. In this frame- 
work, Chandrasekhar [28| observes that the addition 
of a magnetic field onto a flow heated from below 
imparts to the liquid an effective kinematic viscosity 



v e fj oc (H z d) 2 / n m . Only the component of the field 
parallel to the gravity vector is effective. This makes 
the critical Reynolds number for stability proportional 
to 1 + Q/Qo- Using the relation between the torque 
and the critical Reynolds number in the linearly un- 
stable regime (eqs. (|3"5"j) and (|3"7j)). this leads to the 
scaling ([47)) . 



F. Turbulent viscosity 

The turbulent viscosity in the direction perpendic- 
ular to the shear v t can be estimated via the mean 
torque T applied by the fluid to the rotating cylinders 
and the mean turbulent velocity profile. Indeed, this 
torque induces a stress equal to: 



T pv 2 u e 
—r = - — = pvtrdr — , 



(48) 



where A is the area of a cylindrical fluid element at 
radius r, p is the fluid density, and ug is the mean 
azimuthal viscosity. Since a similar formula applies in 
the laminar case, with vt — v 1 one simply gets: 



Vi 



Slam G 
S Gl a „ 



(49) 



Using the expression of Re, Rc and Gi ar . 
thus get the simple expression: 



Sf 2 = /3Sr 



-.2 



(50) 

Here, we have adopted the notation of Richard and 
Zahn 5] to express the turbulent viscosity in unit of 
the typical shear and radius of the flow as j3. This non- 
dimensional parameter encompasses all the interesting 
variation of the turbulent viscosity as a function of the 
radial position r and the control parameters i?e, Rc 
(or rj) and Rn- The radial variation is given through 
the ratio Si am /S as illustrated in section HVBI This 
ratio is one near the boundary and may increase in 
the core of the flow, due to the turbulent shear re- 
duction. All the variation with Rq is through the 
function h which has been empirically determined in 
Section IV Dl and plotted in Fig. [T5] All the variation 
with Re is through Gi(Re,rj) / Re 2 , which can be de- 
termined through torque measurements f Sections IV Al 
and IV Bl) , with a theoretical expression provided in 
Section lVDI for smooth boundaries, and lVEI for rough 
boundaries. The dependence on the curvature is sub- 
tler since it appears in all the above dependencies. 

An example of variation of the dimensionless turbu- 
lent viscosity (3 for -q = 0.72 is provided in Fig. I2TH for 
smooth and rough boundary conditions. One sees that 
at large enough Reynolds number, this function be- 
comes independent of the Reynolds number for rough 
boundary conditions, while it decreases steadily in the 
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smooth boundary cases, due to logarithmic correc- 
tions. This weak Reynolds number variation is in con- 
trast with standard turbulent viscosity prescription, 
based on dimensional consideration a la Kolmogorov. 
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FIG. 20: Influence of Reynolds number on turbulent vis- 
cosity. Case with two rough boundaries at r\ = 0.724, V, 
data from Van den Berg et al.. The continuous line is 
drawn using formula (143[) . Case with two smooth bound- 
aries at r\ = 0.724, A, data from hewis and Swinney. The 
dotted and the dashed-dot lines are drawn using formulae 

(S3- 

Finally, let us compare our results with previ- 
ous results for the turbulent viscosity in rotating 
flows. Using a turbulent closure model of turbulence, 
Dubrulle [52] derived v t =2x 1CT 3 R^ 2 Sf 2 . This for- 
mula reflects the correct behavior in term of Rq (see 
Section l"V CP but fails to reproduce the Reynolds de- 
pendence in the case of smooth boundary conditions. 
For rough boundary conditions, our formula predicts a 
turbulent viscosity going like = 8 X 10~ 3 RQ 2 Sf 2 for 
Ru < -0.5 and c t =2x W^R^SP 2 for R n > 0.5 x 
10 , in the wide gap limit. The formula of Dubrulle 
is therefore in between these two predictions. Richard 
and Zahn @ used Taylor measurements [l6[ to derive 
the value /3 = 1.5 ± 0.5 x 10~ 5 . These measurements 
are performed for Re ~ 5 x 10 , with inner cylinder at 
rest. At r\ — 0.7, one has Rn = 0.4 and from Fig. Q1)J 
h(QA, 0.7) = 0.05. For Si am /S we adopt a value equal 
to 2, as suggested by Section HVB I Finally, from Fig. 
E2 we get R A c G l {Re,Q.7)/Re 2 = 5 x 10~ 4 so that we 
finally obtain j3 = 8 x 10~ 6 , an estimate close to the 
one proposed by Richard and Zahn. 



for the Taylor- Couette flow. This prediction clearly 
indicates the dependencies on the Reynolds number 
and the rotation number. The curvature effect is much 
trickier to isolate at least with the available data, since 
it appears in all the terms of the prescription. Espe- 
cially, on the cyclonic side, where the Rayleigh crite- 
rion depends on the curvature, it is impossible without 
any phenomenological arguments to isolate the cur- 
vature effect from the rotation one within the set of 
data used here. Since we wanted to remain as close as 
possible to the existing data, we decided not to repro- 
duce any phenomenological arguments in the present 
paper (for such an analysis see e.g. Longaretti and 
Dauchot [Til). 

The introduction of new control parameters, which 
rely on the dynamical properties of the flow, rather 
than on its geometry allows us to envision some ap- 
plication of our result to rotating shear flows in general 
even if one should remain cautious with the details of 
the boundary conditions. These new control param- 
eters have a rather general ground, but they remain 
global quantities. It would be interesting to further 
develop this approach, by introducing local dynami- 
cal control parameters, so that in spatially developing 
flows, one could conduct a local study of the stability 
properties. 

In order to validate the above prescription, it would 
definitely be necessary to confront it to more exper- 
imental data. On the anticyclonic side, in the sub- 
critical regime, only one value of the curvature has 
been investigated, so that we have very little idea of 
its influence. In the supercritical regime, an impor- 
tant hypothesis made here was to introduce R s g up and 
to relate it to Rt and R c . Also, we have proposed 
to relate the torque measurements (a difficult exper- 
imental task) to the threshold determination. These 
conjectures should be checked against more data. Fi- 
nally, we have tried to provide some indications on 
the influence of external effects such as stratification, 
or magnetic fields. Clearly the lack of experimental 
data here is such that very little could be done and a 
definitive effort should be conducted in this direction. 

Still, it is to our knowledge the first time that using 
most of the existing experimental studies a practical 
prescription for the turbulent viscosity is proposed. It 
can certainly be improved, but we believe that, even 
at the present level, it can already bring much insight 
into the understanding of some astrophysical and geo- 
physical flows. 



VI. CONCLUSION AND HINTS FOR 
FUTURE WORK 



Acknowledgments 



The present work provides us with a prescription for 
the turbulent viscosity hence the turbulent transport 



We would like to thank L. Marie, F. Hersant and F. 
Busse for fruitful discussions. 



19 



[1] F. Hersant, B. Dubrulle and J-M. Hure, "Turbulence 
in circumstellar disks", A& A accepted (2004). 

[2] D. Richard. Instabilites hydrodynamiques dans les 
ecoulements en rotation differentiate. PhD thesis. 
Universite de Paris VII (2001). 

[3] Y.B. Zeldovich, "On the friction of fluids between ro- 
tating cylinders" , Proc. Roy. Soc. London A, 374, 299 
(1981). 

[4] B. Dubrulle, "Differential rotation as a source of angu- 
lar momentum transfer in the Solar Nebula", Icarus, 
106, 59 (1993). 

[5] D. Richard and J. -P. Zahn, "Turbulence in differen- 
tially rotating flows", Astron. Astrophys., 347, 734 
(1999). 

[6] P.-Y. Longaretti, "On the phenomenology of hydro- 
dynamic shear turbulence", Astrophys. J., 576, 587 
(2002). 

[7] S. K. Bahl, "Stability of viscous flow between two 
concentric rotating porous cylinders", Def. Sci. J. , 
20, 89 (1970). 

[8] K. Min and R. Lueptow, "Hydrodynamic stability 

of viscous flow between porous cylinders with radial 

flow", Phys. Fluids, 6, 144 (1994). 
[9] G. Wendt, "Turbulente Stromung Zwischen Zwei 

Rotierenden Konaxialen Zylindern", Ing. Arch., 4, 

577 (1933). 

[10] D. Coles and C. van Atta, "Measured distortion of a 
laminar circular Couette flow by end effects" , J. Fluid 
Mech. 25, 513, (1966). 

[11] P.-Y. Longaretti, O. Dauchot, "Rotation and curva- 
ture effects in subcritical turbulent shear flows" , sub- 
mitted to Physics of Fluids (2004). 

[12] B. Dubrulle and F. Hersant, "Momentum transport 
and torque scaling in Taylor-Couette flow from an 
analogy with turbulent convection" , Eur. Phys. J. B, 
36, 379 (2002). 

[13] A. Prigent, B. Dubrulle, O. Dauchot and I. Mut- 
tabazi, "The Taylor-Couette flow: the hydrody- 
namic twin of Rayleigh-Benard convection" , accepted 
(2004). 

[14] D.D. Joseph, "Stability of Fluid Motions", I, Spriger 
tracts in Natural Philosophy, Vol. 27, Springer Verlag 
(1976). 

[15] O. Dauchot, P. Manneville, "Local vs. Global Con- 
cepts in Hydrodynamic Stability theory", J. PHys. 
II, France, 7 pp. 371-389 (1997). 

[16] G.I. Taylor, "Fluid friction between rotating cylin- 
ders", Proc. Roy. Soc. London, A 157, 546 (1936). 

[17] N. Tillmark and P.H. Alfredsson, "Experiments on ro- 
tating plane Couette flow" , in Advances in Turbulence 
VI, Gavrilakis, L. Machiels and P. A. Monkewitz eds, 
Kluwer Academics Publishers, p391, (1996). 

[18] A. Prigent, G. Gregoire, H. Chate, O. Dauchot and 
W. van Saarloos, "Long wavelength modulation of 
turbulent shear flows" , Phys. Rev. Letters, 89, 014501 
(2002). 

[19] G.S. Lewis and H.L. Swinney, "Velocity structure 
function, scaling and transitions in high-Reynolds 
number Couette) Taylor flow", Phys Rev. E 59 5457 
(1999). 

[20] M.M Couette, "Etudes sur le frottement des liquides" , 



Ann. Chim. Phys., 6, 433 (1890). 
[21] G.I. Taylor, "Stability of a viscous liquid contained 

between two rotating cylinders", Phil. Trans. Roy. 

Soc. London A, 223, 289, (1923). 
[22] R.J. Donnelly and D. Fultz, "Experiments on the sta- 
bility of spiral flow between rotating cylinders" , Proc. 

Nat. Acad. Set., 46, 1150 (1960). 
[23] H.A. Snyder, "Stability of rotating Couette flow. II. 

Comparison with numerical results", Phys. Fluids, 

11, 1599 (1968). 
[24] P. Bradshaw, "The analogy between streamline curva- 
ture and buoyancy in turbulent shear flow", J. Fluid 

Mech., 36, 177 (1969). 
[25] T.J. Pedley, "On the instability of viscous flow in a 

rapidly rotating pipe", J. Fluid Mech., 35, 97 (1969). 
[26] D.K. Lezius and J. P. Johnston, "Roll-cell instabilities 

in rotating laminar and turbulent channel flows", J. 

Fluid Mech., 77, 153 (1976). 
[27] V.A. Romanov, "Stability of plane-parallel Couette 

flow", Functional Anal. Applic, 7, 137 (1973). 
[28] S. Chandrasekhar, "The stability of non-dissipative 

Couette flow in hydromagnetics" , Proc. Nat. Acad. 

Sci. 46, 253 (1960). 
[29] A. Esser and S. Grossmann, "Analytic expression for 

Taylor-Couette stability boundary", Phys. Fluids, 8, 

1814 (1996). 

[30] D.P. Lathrop, J. Fineberg, H.L. Swinney, "Transition 
to shear driven turbulence in Couette- Taylor flow", 
Phys. Rev. A, 46, 6390 (1992). 

[31] O. Dauchot and F. Daviaud, "Finite amplitude per- 
turbation and spot growth mechanism in plane Cou- 
ette flow", Phys. Fluids 7, 335 (1994). 

[32] K. Park, G.L. Crawford, and R.J. Donnelly, "Deter- 
mination of transition in Couette flow in finite geome- 
tries", Phys. Rev. Letters, 47, 1448 (1981). 

[33] E.P. Velikhov, "Stability of an ideally conducting liq- 
uid flowing between cylinders rotating in a magnetic 
field", Sov. Phys. JETP 9, 995, (1959). 

[34] R.J. Donnelly and M. Ozima, "Experiments on the 
stability of flow between rotating cylinders in the pres- 
ence of a magnetic field", Proc. Roy. Soc. A, 266, 272, 
(1962). 

[35] S.A. Balbus and J.F. Hawley, "A powerful local shear 
instability in weakly magnetized disks. I. Linear anal- 
ysis", Astrophys. J., 376, 214 (1991). 

[36] G. Riidiger, M. Schultz and D. Shalybkov, "Linear 
magnetohydrodynamic Taylor-Couette instability for 
liquid sodium", Phys. Rev. E, 67, 046312 (2003). 

[37] A. P. Willis and C.F. Barenghi, "Magnetic instability 
in a sheared azimuthal flow", A&A, 388, 688 (2002). 

[38] L. Howard and A. Gupta, "On the hydrodynamic and 
hydromagnetic stability of swirling flows" J. Fluid 
Mech.,14, 463 (1962). 

[39] Dubrulle, B. and Knobloch, E. "On instabilities in 
magnetized accretion disks", A&A,256, 673 (1992). 

[40] M.J. Molemaker, J.C. McWilliams and I. Yavneh, 
"Instability and equilibration of centrifugally stable 
stratified Taylor-Couette flow" Phy. Rev. Letters, 86, 
5273 (2001). 

[41] B. Dubrulle, L. Marie, Ch. Normand, F. Hersant, 
D. Richard and J-P. Zahn, "An hydrodynamic shear 



instability of stratified disks", submitted to A& A 
(2003). 

[42] I. Yavneh, J.C. McWilliams and M. J. Molemaker, 
"Non-axisymmetric instability of centrifugally stable 
stratified Taylor-Couette flow" J. Fluid Mech., 448, 
1 (2001). 

[43] E.M. Withjack and C.F. Chen, "An experimental 
study of Couette instability of stratified fluids", J. 
Fluid Mech., 66, 725 (1974). 

[44] B.M. Boubnov and E. Hopfinger, "Stability of the 
Couette flow between Two independently rotating 
cylinders in a stratified fluid", Physics-Doklady 42, 
312 (1997) 

[45] J-C Chen and J-Y. Kuo, "The linear stability of 
steady circular Couette flow with a small radial tem- 
perature gradient", Phys. Fluids, A2, 1585 (1990). 

[46] I. Mutabazi, A. Goharzadeh and F. Dumouchel, "The 
circular Couette flow with a radial temperature gradi- 
ent" in 12 th International Couette-Taylor Workshop, 
September 6-8, 2001, Evanston, IL, USA, (2001). 

[47] HA. Snyder and S.F. Karlsson, "Experiment on the 
stability of Couette motion with a radial thermal gra- 
dient", Phys. Fluids 7, 1696 (1964). 

[48] F.H. Busse, "A property of the energy stability limit 
for plane parallel shear flow" , Arch. Rat. Mech. Anals, 
47, 28 (1972). 

[49] F.H. Busse, "Bounds for properties of complex sys- 
tems", in "Non-linear Physics of Complex Systems", 
G. Parisi, S.C. Muller, W. Zimmermann (eds), Lec- 
ture Notes in Physics, 476, 1 (1996). 

[50] F.H. Busse, "Bounds for turbulent shear flows", J. 
Fluid Mech., 41, 219 (1970). 

[51] L. Jacquin, O. Leuchter, C. Cambon and J. Math- 
ieu, J. Fluid Mech., "Homegeneous turbulence in the 
presence of rotation", 220 1 (1990). 

[52] B. Dubrulle and L. Valdettaro, "Consequences of ro- 
tation in energetics of accretion disks", A&A, 263, 
387 (1992). 

[53] T.H. van den Berg, C.R. Doering, D. Lohse and D.P. 
Lathrop, "Smooth and rough boundaries in turbulent 
Taylor-Couette flow", Phys. Rev. E, 036307 (2003). 

[54] O. Cadot, Y. Couder, A. Daerr, S. Douady, A. Tsi- 
nober, "Energy injection in closed turbulent flows; 
stirring through boundary layers versus inertial stir- 
ring", Phys. Rev. E , 56, 427 (1997). 

[55] Werne, J., "Turbulent convection: what rotation has 
taught us?", in Geophysical and Astrophysical Con- 
vection, Eds Fox and Kerr, Gordon Breach, 221, 
(1996). 

[56] B. Dubrulle, "Logarithmic correction to scaling in tur- 
bulent thermal convection", Eur. Phys. J. B, 21, 295 
(2001). 

[57] B. Dubrulle "A turbulent closure model for thin ac- 
cretion disks", A&A,2QG, 592 (1992). 



21 



TABLE II: Notations 



Tj ■■ 
h 



Vijl 



Superscript and subscript conventions 

any flow variable (e.g., component of velocity) 
laminar part of X \ ~ 2 / 

mean part of X £ J ~ . 

fluctuating part of X R±-±Vd/v 
typical value of X 
relates to inviscid flows 
relates to Plane Couette flows 
relates to Taylor-Couette flows 

relates to subcritical flows . . 

relates to supercritical flows [L\ = 

relates to cyclonic flows [T] = 

relates to anti-cyclonic flows R e = 

inertial frame cartesian components of X, 
inertial frame cylindrical 3 components o] 
rotating frame cylindrical" components ®fiX. 

Hydrodynamical quantities 
position vector 
inertial frame velocity vector 
rotating frame velocity vector 
fluid pressure, generalized pressure 
angular velocity 

angular velocity of the rotating frame 
velocity shear 

du l ™ I 'dy; du x / 'dy (plane Couette flow) 



radius ratio (dimensionless measure of the gap) 

cylinders height 

Dimensionless quantities 

Standard adimensionalization 
unit of length 

unit of time 

Reynolds number of moving plates (plane Couette flo' 
Reynolds number (plane Couette flow) 
Reynolds number of rotating cylinders (Taylor- Couett 
Reynolds number (Taylor-Couette flow) 
Dynamical adimensionalization 
(Taylor-Couette flow) 
unit of length 
unit of time 
Reynolds number 
rotation number 

rotation number at marginal stability (Re = co) 
-curvature number 

Local dynamical ratios 
(Taylor-Couette flow) 
Local "Reynolds" ratio 
Local "rotation" ratio 
Local "curvature" ratio 

Transition Reynolds numbers 
first supercritical linear transition 
minimal Reynolds number for self-sustained turbulene 
Transition to "hard" turbulence (as traced by torques 

Symbol conventions in graphs 
cyclonic data 
anti-cyclonic data 
Wcndt (1933) data 

Taylor (1936) data 

Continued on next page 



0, ♦ Tillmark and Alfredsson (1996) data 

o, • Richard (2001) data 

V, T Lewis and Swinney (1999) data 

a The corrcspondancc between cartesian and cylindrical axes is (x —0, — 0), and 



X 

X 
X' 
X 

x°° 

X vCi 

x TCf 

-y-sub 
yrsup 

X+ 

x~ 

X — (X X , Xy , X- 

X = (X r , Xg , X z 

X — (X r , X<t,, X; 



Re = R+ - R- 
Re = \R -RA 



d 

s- 1 

Sd 2 /v _ 
2Q. r j I S 



x 
u 
w 

P ,7T 

fi = Ug/r 

firf 
s 

L = (as x u) z 
T 

G = T/(phu 2 - 



r„ = 

r c = 



ad vection — shear 
viscous clissipatior 

Coriolis 
acl vection — shear 
curvature 



advecti* 



-shear 



rdn lam /dr;rdn/dr (Taylor-Couette flo^ 
= r fl specific angular momentum 
torque 

adimensionalized torque 



□ 



(Geometric and physical quantitiesU, 



0,o, A 



v 
v t 
P 

±V 



d = r a — n 



General 
kinematic viscosity 
turbulent viscosity 
mass density 

Couette flow 
bounding plate velocity 
gap between bounding plates 
Taylor-Couette flow 
Inner, outer cylinder radii 
Inner, outer cylinder angular velocity 
gap 
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